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Resonant magnetic x-ray diffraction experiments were carried out on the stacked triangular 
lattice antiferromagnet GdPd2Al3. The experiments revealed an expected initial collinear c- 
axis order at Xni followed by an additional in-plane order at Tn2, while at the same time 
we found that the ground state is a helically ordered state of a very long incommensurate 
period of approximately 700 A. The distribution of K-domains was highly anisotropic, and 
the domain with the modulation vector normal to the surface of the crystal was ascendant. 
Low-field magnetization is discussed on the basis of the observed incommensurate magnetic 
structure. 
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1. Introduction 

Geometrical frustration has occupied a central position in condensed matter physics for 
decades. 1 Stacked triangular lattice antiferromagnets were extensively investigated at the 
early stages of the history, particularly on the quasi-one-dimensional hexagonal ABX% anti- 
ferromagnets. Although these magnets do not posses liquid-like ground states, 2-4 which are a 
hallmark of highly frustrated magnets, and undergo phase transitions to long-range ordered 
states, the novel ordered states and the distinctive nature of phase transitions found in these 
magnets illuminate prominent aspects of various frustration-related phenomena. 5-7 

Phase transitions of Heisenberg spins with weak Ising anisotropy on a triangular lattice 
are very well investigated. 8 The characteristic features of these magnets are summarized in 
the following three points, (i) At zero field, successive phase transitions of the z- and xy- 
components of the spins occur. As the temperature decreases, the z-component first enters 
the long-range ordered state at Tni as shown in Fig. 1(a), and then the xy-components exhibit 
long-range ordering at Tn2- Finally, a slightly distorted 120° structure (Fig. 1(b)) is completed 
below Tn2- (ii) Since the canting angle a is smaller than 60°, the total sum of the magnetic 
moments in a triangle does not cancel out. Accordingly, a net moment along the z-axis is 
observed in a single triangular plane, (iii) When the magnetic field is applied along the z-axis, 
the magnetization is maintained at a constant value at a certain range of magnetic fields. The 
value is exactly one-third of the saturation magnetization, and the phenomenon is known as 
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Fig. 1. Magnetic structures of Heisenberg triangular lattice antifcrromagnct with weak Ising 
anisotropy. Projection onto the (110) plane, (a) Between Tni and Tn2 and (b) below Tn2- The 
canting angle a is 60° for the pure Heisenberg case and decreases as the weak Ising anisotropy 
increases, (c) Example of incommensurate structures. 

the "one-third (magnetization) plateau". 

The first point, successive phase transitions, was experimentally confirmed by Clark and 
Moulton more than 30 years ago in CsNiCi3, 9 which is one of the above-mentioned ABX3 
compounds. However, the predominant antiferromagnetic coupling along the chain of these 
compounds inhibits the appearance of the second and third points, namely a small net moment 
along the z-axis and the one-third magnetization plateau. The experimental observation of the 
latter two points was therefore delayed until the recent discovery of the prototypical compound 



form a stacked triangular lattice. The magnetic properties of this compound were investi- 
gated in detail by Kitazawa and coworkers. 10 ' 12 The magnetic susceptibility well follows the 
Curie- Weiss law, and the Curie temperature is reported to be —30 K. Hence, the dominant 
interactions are antiferromagnetic. The specific heat measurements demonstrated successive 
phase transitions at Tni = 16.8 K and Tn2 = 13.3 K. The magnetization at low fields indicated 
a weak ferromagnetic moment along the c-axis. Therefore, it is expected that triangular layers 
with small net moments couple ferromagnetically with each other. At high fields, the one-third 
magnetization plateau was observed between 6.2 T and 11.8 T, when only the magnetic field 
parallel to the c-axis was applied. All of these pieces of evidence are clearly in line with the 
criteria of a Heisenberg stacked triangular lattice antiferromagnet with weak Ising anisotropy. 
Nevertheless, the magnetic structures of GdPd2Ak> have not yet been investigated, since neu- 
tron diffraction experiments of Gd compounds are fairly difficult to perform owing to the large 
absorption cross section of natural Gd for neutrons. No microscopic measurements have been 
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reported thus far, with the exception of a Mdssbauer spectroscopy experiment, which suggests 
a non-collinear magnetic structure of more than two sublattices. 11 The order parameters at 
Tni and Tn2 have not yet been revealed. 

In this paper, we carried out resonant magnetic x-ray diffraction experiments on 
GdPd2Al3. The temperature dependence of rotated and unrotated signals with respect to 
the incident x-ray polarization unambiguously illustrated successive phase transitions of the 
c-axis and a6-plane components at Tni and ?N2, respectively. An unexpected outcome was 
that the magnetic structure below Tn2 was incommensurately modulated. The period was 
found to be about 130 times as long as the lattice constant a. In addition, it was found that 
three domains with respect to the direction of the modulation wave vector (K-domains) were 
distributed extremely unequally. We also discuss the low-field magnetization from the view 
point of this long-period incommensurate magnetic structure. 

2. Experimental 

Resonant x-ray diffraction experiments were carried out at beamline BL22XU in SPring- 
8. The photon energy was tuned near the L2 absorption edge of Gd (7.930 keV). A single 
crystal of GdPd2Ah> was grown by the Czochralsky pulling method in a tetra-arc furnace. 
The sample was cut into a parallelepiped of 4 x4 x2 mm 3 in volume, and a (110) surface was 
polished. The sample was attached to the cold head of a conventional closed-cycle refrigerator, 
which was mounted on a conventional four-circle diffractometer with a horizontal scattering 
plane. We set the c-axis of the sample perpendicular to the scattering plane. The orientation 
of the sample was determined using the 110 and 111 reflections, and we confirmed that the 
angle between the c-axis and the normal vector of the scattering plane was less than 7°. The 
polarization of the scattered x-rays was analyzed with respect to whether it was normal (a) 
or parallel (tt) to the scattering plane using the 006 reflection of a pyrolytic graphite crystal. 
Since the incident polarization was parallel to the scattering plane (ir), the scattered x-rays 
were separated into rotated (tt to a') and unrotated (tt to tt') channels, where the prime sign 
represents scattered x-rays. In order to prevent the beam from heating the sample surface, 
we reduced the incident photon flux by a factor of 7. The experiments were carried out twice. 
In the first experiment, the mosaic width of the 110 reflection of the sample was about 0.07° 
full width at half maximum (FWHM) , indicating the high quality of the crystal. Prior to the 
second experiment, the sample surface was carefully polished again, and the mosaic width was 
improved to 0.015° FWHM. The data provided in the rest of this paper were obtained from 
the repolished sample, unless explicitly specified otherwise. 

3. Results 

First, we show the peak intensity of the magnetic reflection ||0 as a function of pho- 
ton energy at 4 K, as well as the fluorescence spectrum, in Fig. 2. Absorption correction 
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Fig. 2. (color online) Energy dependence of the peak intensity of a magnetic Bragg reflection at 
(|, g,0). The fluorescence spectrum is also shown. Absorption correction is not performed. The 
diffraction intensity resonantly increases near the Gd L2 absorption edge 7.931 keV. 



was not performed. The fluorescence spectrum indicates that the Gd L2 absorption edge is 
7.931 keV. The diffraction intensity was very weak far from the L2 absorption edge, while an 
enormous resonant enhancement of the intensity was observed near the main edge, implying 
that the resonance is ascribed to electric dipole (El) transitions (2p — > 5d). Subsequent mea- 
surements were performed at the peak energy (7.932 keV). We confirmed that the modulation 
wave vector q M is only (g, g, 0) in the first experiment. Reciprocal lattice scans from (|, |, 0) 
to (|,|,1) and from (0.4,0.4,0) to (1.7,1.7,0) yielded four magnetic reflections at (|,|,0), 
(|,|,1), (3,3,0), and (|,|,0). No other modulation wave vectors were detected. We also 
surveyed an entire Brillouin zone by using an area detector and confirmed the above result. 13 

The main result was obtained from the temperature dependence of the ||0 reflection. In 
Fig. 3(a), the diffraction intensity integrated along the radial (9-29) direction is shown as a 
function of temperature for both the tt-tt' and tt-g' channels. The data unambiguously illus- 
trate that the tt-tt' channel appears below 17.3 K (Tni), whereas the tt-o' channel is observed 
only below 14.3 K (Tsf2)- The slight discrepancies (about 1 K) between the observed transition 
temperatures and the values cited in the literature are probably due to radiation heating of 
the thermometer. In the intermediate phase between Tni and Tsf2, only the scattering process 
from 7T to tt' is allowed. 

Here we briefly describe the resonant x-ray scattering amplitude / rcs of Gd ions. According 
to ref. 14, / rea is proportional to 

/ res oc C ef ■ e ; + iCi{e\ x e ; ) ■ m + C 2 e\ O e i5 (1) 

where Co, C±, and C2 are energy-dependent constants, m is the magnetic moment, and £j 
and ef are the polarization vectors of the incident and scattered x-rays, respectively. The 
symmetric second-rank tensor O describes the anisotropy of the Gd 5d orbital caused by 
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Fig. 3. (color online) (a) Integrated intensity of the rotated (tt-ct') and unrotated (tt-tv') channels of 
the ||0 magnetic reflection as a function of temperature. The rotated and unrotated channels 
sense the magnetic moments normal and parallel to the c-axis, respectively. See the text for details, 
(b) Deviations of the peak position from the commensurate position (A) for the 110, ||0, and 
||0 reflections. The unit is a reduced lattice unit. The unrotated channels are depicted for the 
magnetic reflections. The deviations of the magnetic reflections clearly develop below Tn2- 

the anisotropic crystal environment. We neglect a small magnetic contribution in O, which 
is proportional to m 2 . 14,15 The first term is an ordinary anomalous scattering factor. The 
second term corresponds to magnetic scattering, and the last term causes the anisotropy of 
the tensor of susceptibility (ATS) scattering. Since all Gd ions are equivalent in GdPd2Al3, 
the first term does not contribute to super-lattice reflections. It is also expected that the last 
term, ATS scattering, does not produce any intensity at super-lattice positions. This is due 
to the fact that Gd ions have no quadrapole moments. Hence, it is reasonable to consider 
that the magnetic ordering does not introduce anisotropic lattice distortions around the Gd 
ions. Therefore, the second term, magnetic scattering, is the only term to be considered in 
GdPd 2 Al 3 . 

As already mentioned, it polarization is parallel to the scattering plane. In addition, we set 
the c-axis of the crystal perpendicular to the scattering plane. Since the resonant x-ray mag- 
netic scattering amplitude includes an outer product of e; and ef, the tt-tt' channel observes 
magnetic moments perpendicular to the scattering plane, which are the c-axis component of 
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Fig. 4. (color online) Intensity map on the (H,K,0) plane for the ||0 reflection at 5.8 K. The intensity 
is plotted on a logarithmic scale, and the units are counts per second. Vertical broken lines indicate 
the positions of (H + K)/2 = 4/3 - 36, 4/3 - S, 4/3, and 4/3 + 26, from left to right. A very weak 
peak might exist at the lower right of the strong peak. Higher harmonic satellites (at 2qM and 
3<?m) were not observed. 

the magnetic moments. In contrast, the scattering amplitude of the ir-a' channel is propor- 
tional to magnetic moments parallel to the incident x-rays, which are the magnetic moments 
in the afr-plane. By applying this characteristic polarization dependence of resonant magnetic 
x-ray diffraction, the temperature dependence of the magnetic reflection shown in Fig. 3(a) 
directly leads to the conclusion that the phase transition at Tni is an ordering of the z- 
component of the magnetic moments and that the xy-components are paramagnetic until the 
temperature drops below Tn2, in good accordance with the expected behavior of Heisenberg 
triangular lattice antiferromagnets with weak Ising anisotropy. 

A minute exploration of the data, however, detects unexpected behavior as well. Firstly, 
we found that the transition at Tn2 is a commensurate-to-incommensurate transition. The 
peak positions of the ||0 and ||0 magnetic reflections moved slightly along the radial (9-26) 
direction below Tn2- We show the deviations of the peak position from the commensurate 
position (A) for the tt-7t' channels of two magnetic reflections as a function of temperature in 
Fig. 3(b), where we denote the wave vectors of the magnetic reflections as q = (h+A, h + A, 0) 
and h = | or |. We also show A for the 110 lattice reflection as a reference. The commensurate 
positions were estimated at 15 K. The deviations for the ||0 and ||0 reflections developed 
positively and negatively below Tn2 like an order parameter of a second-order phase transition, 
respectively, whereas the deviation for the 110 reflection remained zero throughout the entire 
range of temperature in the measurements. The magnetic modulation wave vector q M is 
thus expressed as (| — S, | — 6,0). 6 was evaluated to be 0.0013 at 5 K. About 130 lattice 
points are therefore included in the period which is approximately 700 A. An example of 
incommensurate magnetic structures is shown in Fig. 1(c). 
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Secondly, the distribution of domains was not uniform. In the incommensurate phase, there 
were three domains with respect to the direction of the modulation vectors (the so-called K- 
domains). For instance, the g-g-0 reflection splits into three Bragg points at (| — <5, g — 5,0), 
(1-^1+ 25,0), and (| + 25, | — 5, 0). We measured the intensity map on the (H,K,0) 
reciprocal-lattice plane around the g|0 reflection at 5.8 K. The result is depicted in Fig. 4. 
There should be three domains. However, only one domain is visible. Although a very weak 
domain might exist at (g — 5, g + 25, 0), most of the intensity concentrates at the primary 
peak. The modulation vector of the strongest domain was normal to the crystal surface. We 
believe that residual strains introduced by the polishing of the (110) surface select the observed 
domain. In the first experiment, the domain distribution was also quite anisotropic. However, 
we observed another domain rather clearly, and hence we consider that repolishing reduced 
the volume of minor domains. Figure 4 also illustrates that there are no higher harmonic 
satellites. The absence of a peak at 2qm suggests that the m 2 term in eq. (1) is actually 
negligible in GdPd2Ai3. 15 No observable intensity at 3qM corresponds to a small squaring up 
of the modulation, indicating that the magnetic structure is described only by the fundamental 
wave vector qM- 

The residual strains might affect the direction of the magnetic moments. Domains with 
respect to the direction of magnetic moments are referred to as S-domains. In GdPd2Ah>, the 
magnetic moments share a single a*c-plane in an S-domain. Since three equivalent a*c-planes 
exist in a hexagonal lattice, there are three S-domains. The distribution of S-domains can be 
estimated from the ratio of the intensity of the tt-tt' channel to that of the ir-a' channel if the 
canting angle a (see Fig. 1(b)) is known, where we assume the commensurate structure for 
the sake of simplicity, a is evaluated to be 59° from the weak ferromagnetic moment along 
the c-axis at T. Although this value requires a slight correction, as seen in the discussion 
section, it is certain that a is very close to 60° and hence we assume a=60°. The ratios of 
tt-g' to tt-tt' reached constant values 1.36 and 5.1 at low temperatures for the ||0 and ||0 
reflections, respectively. This result leads to the conclusion that about 70% of domains have a 
spin plane normal to the modulation vector. In contrast, the equal distribution of S-domains 
provided a reasonable account in the first experiment. 

4. Discussion 

4-1 Magnetic structure 

Throughout these experiments, we found that there is only one magnetic wave vector 
q M = (| — 5, | — 5, 0) in GdPd2Ai3. In addition, we observed the magnetic moments parallel 
to the c-axis My and the magnetic moments perpendicular to the c-axis M±. Therefore the 
magnetic moment M below Tn2 at the position r is given by 

M(r) = Mm sin(q M r + fo) + M ± sm(q M r + (f> + <f>). 
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Fig. 5. Low-field magnetization at 5 K when magnetic field parallel to the c-axis is applied. (Data 
taken from ref. 10.) 

At <f> = it/2, the magnetic structure is a helically modulated structure such as that shown in 
Fig. 1(c). On the other hand, when 4> = 0, the magnetic structure is a sinusoidally modulated 
one, which inclines from the c-axis by 6 C = tan _1 (|M"_i_|/|M"|| |). Although it is not easy to 
determine the value of cf> experimentally, 16 it can be reasonably inferred that (f> is ir/2. Firstly, 
a sinusoidally modulated structure is unfavorable with regard to entropy at low temperatures. 
Fully polarized magnetic moments are expected for such localized magnets. Secondly, an 
inclined structure is not stable. The origin of the anisotropy must be a dipole-dipole interaction 
since the Gd ion has no orbital moment. Owing to the ferromagnetic coupling between the 
layers along the c-axis, the dipole-dipole interaction behaves as effective Ising anisotropy. If a 
sinusoidally modulated structure is realized, 6 C tends to be zero. Accordingly, we conclude that 
the magnetic structure below Tn2 is a helically and incommensurately modulated structure. 
4-2 Incommensurate modulation 

Normally, an incommensurate structure is ascribed to competition between the nearest- 
neighbor interaction and further- neighbor interactions. 17 Long-range RKKY interactions are 
relevant for the metallic GdPd2Al3. Hence, further- neighbor interactions can be comparable 
in magnitude to the nearest-neighbor interaction. It is reasonable to assume that the observed 
incommensurate structure is stabilized by second- and third-neighbor antiferromagnetic inter- 
actions. These isotropic interactions, however, do not explain why the intermediate phase is 
commensurate. The concomitant development of the incommensurability 8 and the in-plane 
component Mj_ below Tn2 shown in Fig. 3 strongly suggests that M± is responsible for the 
incommensurate structure through anisotropic interactions. 

In-plane incommensurate magnetic structures in stacked (hexagonal) triangular lattice 
antiferromagnets were observed in RbFeCl3. 18 The phase transitions in RbFeCl3 are quite 
peculiar. There are two (linearly-polarized) incommensurate structures at high temperatures, 
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and the low-temperature structure is a commensurate coplanar one. These incommensurate 
magnetic structures and the novel ordering process are well interpreted by introducing dipole- 
dipole interactions between the in-plane magnetic moments. 19 The exchange energy has a min- 
imum at the K point (= (|^, |^ , 0)), while the dipole-dipole interactions between the in-plane 
magnetic moments prefer other ordering wave vectors. Hence, stable points appear near the K 
point, and incommensurate structures are realized. In contrast, dipole-dipole interactions be- 
tween the magnetic moments parallel to the c-axis are stable at the K point. Therefore, dipole- 
dipole interactions provide a reasonable account of incommensurate-commensurate transitions 
at In2- Detailed calculations will be published elsewhere. 
4-3 Low-field magnetization 

In the rest of the paper, we discuss the low- field magnetization of GdPd2Al3. As shown in 
Fig. 5, it appears that the magnetization along the c-axis has a small ferromagnetic moment. 
The extrapolation of the magnetization curve above 0.02 T intersects about 0.22 ^e/Gd at 
T. If the magnetic structure is commensurate, this weak net moment is in good agreement 
with the expectation that a distorted 120° structure due to the weak Ising anisotropy gives rise 
to a small ferromagnetic moment in a single triangular plane since the triangular planes couple 
ferromagnetically with each other, as experimentally confirmed by the zero z-component of 
the modulation wave vector. However, the observed magnetic structure is incommensurate. 
Since the magnetic structure is very close to the 120° structure, a local triangle of spins may 
have net moment, for instance, along the +z direction. However, at a position a half-period 
away, the net moment of the triangle points to the —z direction. Therefore, the triangular 
plane has no ferromagnetic moment. This discrepancy can be resolved by inserting domain 
walls. If the spins rotate clockwise in the first half-period of the incommensurate structure and 
rotate counterclockwise in the next half-period, the total magnetic moments in all half-period 
domains point in the same direction. It is likely that this domain structure is stable in the 
presence of magnetic fields. In fact, this domain structure can be regarded as a fan structure, 
which was introduced by Nagamiya et al. 20 in the helical-to-fan transition of a helimagnet 
under magnetic fields. The transition field to the fan structure is considered to be low, since 
the density of the domain walls is extremely low owing to the long incommensurate period. As 
the magnetic field increases, the amplitude of the incommensurate structure decreases, and 
eventually the commensurate structure is achieved. 

5. Conclusions 

We investigated the magnetic structures and phase transitions of the Heisenberg triangular 
lattice antiferromagnet GdPd2Al3 by means of resonant magnetic x-ray diffraction. Utilizing 
the characteristic polarization dependence of resonant magnetic x-ray diffraction, we revealed 
that only the c-axis component of the magnetic moments exhibits long-range ordering between 
Tni and Tn2 and that the helical structure is completed by the ordering of the a6-plane 
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components below Tn2- The magnetic structure below Tn2 was found to be incommensurate 
and is most likely stabilized by the dipole-dipole interactions between the in-plane components 
of the magnetic moments. We also found that both the K-domains and the S-domains were 
not evenly distributed. It is likely that the residual strains caused by polishing the surface 
produce these asymmetric distributions. The low-field magnetizations were discussed from 
the viewpoint of the observed incommensurate magnetic structures, and we infer that the fan 
structure is stabilized in the magnetic field. 
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